Constraints on massive sterile plus active neutrino species in non minimal cosmologies 
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Cosmological measurements are affected by the energy density of both active and sterile massive 
neutrinos. We extend here a recent analysis of current cosmological data to non minimal cosmologies. 
Several possible scenarios are examined: a constant w 7^ — 1 dark energy equation of state, a non 
flat universe, a time varying dark energy component and coupled dark matter dark energy universes 
or modified gravity scenarios. When considering cosmological data only, (3+2) massive neutrino 
models with ~ 0.5 eV sterile species are allowed at 95% CL. This scenario has been shown to 
reconcile reactor, LSND and MiniBooNE positive signals with null results from other searches. Big 
Bang Nucleosynthesis bounds could compromise the viability of (3+2) models if the two sterile 
species are fully thermalized states at decoupling. 

PACS numbers: 98.80.-k 95.85.Sz, 98.70.Vc, 98.80.Cq 



I. INTRODUCTION 

Solar, atmospheric, reactor, and accelerator neutrinos 
have provided compelling evidence for the existence of 
neutrino oscillations, implying non-zero neutrino masses 
(see Ref. [l| and references therein). The present data 
require the number of massive neutrinos to be equal or 
larger than two, since there are at least two mass squared 
differences (Am 2 tmos and Am 2 olar ) driving the atmo- 
spheric and solar neutrino oscillations respectively. Un- 
fortunately, oscillation experiments only provide bounds 
on the neutrino mass squared differences, i.e. they are 
not sensitive to the overall neutrino mass scale. 

Cosmology provides one of the means to tackle the ab- 
solute scale of neutrino masses. Neutrinos can leave key 
signatures in several cosmological data sets. The amount 
of primordial relativistic neutrinos changes the epoch of 
the matter- radiation equality, leaving an imprint on 
both Cosmic Microwave Background (CMB) anisotropics 
(through the so-called Integrated Sachs- Wolfe effect) and 
on structure formation, while non relativistic neutrinos in 
the recent Universe suppress the growth of matter density 
fluctuations and galaxy clustering, see Ref. 0. Cosmol- 
ogy can therefore weigh neutrinos, providing an upper 
bound on the sum of the three active neutrino masses, 
Y,m v ~ 0.58 eV at 95% CL [|. The former bound is 
found when CMB measurements from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) are combined with 
measurements of the distribution of galaxies (SDSS-II- 
BAO) and of the Hubble constant H (HST) 1 in the 
assumption of a flat universe with a cosmological con- 
stant, i.e. a ACDM cosmology. 

There is no fundamental symmetry in nature forcing a 
definite number of right-handed (sterile) neutrino species, 
as those are allowed in the Standard Model fermion con- 
tent. Indeed, cosmological probes have been extensively 



1 For other recent analyses, see also Refs. 0, 



used to set bounds on the relativistic energy density of 
the universe in terms of the effective number of neutrinos 
Nf (see, for instance, Refs. 0-i]. If the effective number 
of neutrinos Nf is larger than the Standard Model pre- 
diction of Nf — 3.046 at the Big Bang Nucleosynthesis 
(BBN) era, the relativistic degrees of freedom, and, con- 
sequently, the Hubble expansion rate will also be larger 
causing weak interactions to become uneffective earlier. 
This will lead to a larger neutron-to-proton ratio and will 
change the standard BBN predictions for light element 
abundances. Combining Deuterium and 4 He data, the 
authors of Ref. [6j found Nf = 3.lt\i at 95% CL. 

Models with one additional ~ 1 eV massive sterile 
neutrino, i.e. the so called (3+1) models, were intro- 
duced to explain LSND short baseline (SBL) antineu- 
trino data [lOj by means of neutrino oscillations. A much 
better fit to SBL appearance data and, to a lesser ex- 
tent, to disappearance data, is provided by models with 
two sterile neutrinos (3+2) [ll[ [l2[ which can also ex- 
plain both the MiniBooNE neutrino 1131 and antineutrino 
data [HI if CP violation is allowed |l5j. More recently, 
a combined analysis including the new reactor antineu- 
trino fluxes [H|,[l3] has shown that (3+2) models provide 
a very good fit to short baseline data [l8|. While these 
models with extra sterile species show some tension with 
BBN bounds on Nf , the extra sterile neutrinos do not 
necessarily have to feature thermal abundances at de- 
coupling, see Refs. [l9|,[2(|, where the usual full thermal- 
ization scenario for the sterile neutrino species was not 
assumed. Up to date cosmological constraints on massive 
sterile and active neutrino species have been presented in 
Refs. [HI, HI] in the context of a ACDM universe. It is 
well known that bounds on active neutrino species are 
relaxed if the dark energy equation of state is different 
from —1 [23M25I ] and/or interactions between the dark 
matter and dark energy sectors are switched on [2(| [53] ■ 
In the same line, the authors of Ref. [HI have found that 
in models with non zero curvature and two extra sterile 
neutrinos the cosmological constant scenario is ruled out 
at95%CL. 
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Here we extend the minimal cosmological scenario 
considered in our previous study [22[ and compute the 
bounds on the masses of the active and the sterile neu- 
trino states as well as on the number of sterile states 
in the presence of a constant equation of state w ^ 1, 
a time varying dark energy fluid, a non vanishing curva- 
ture component and interactions among the dark sectors. 
The paper is organized as follows. Section [TT] describes 
the details of the analysis carried out here, including the 
cosmological parameters and datasets. The four different 
cosmological scenarios explored here are analyzed and the 
most important degeneracies among the neutrino param- 
eters are carefully explored. Section IlIII summarizes our 
main results and conclusions. 



II. COSMOLOGICAL CONSTRAINTS 

Here we present the constraints from current data on 
the active neutrino masses and on the sterile neutrino 



thermal abundance and masses in different cosmologi- 
cal scenarios. We have modified the Boltzmann CAMB 
code p^ l incorporating the extra massive sterile neu- 
trino parameters and extracted cosmological parameters 
from current data using a Monte Carlo Markov Chain 
(MCMC) analysis based on the publicly available MCMC 
package cosmomc (30j . We consider here four possible 
scenarios: the wCDM model in which we include the 
possibility of a dark energy equation of state parameter 
w different from -1; the w(o)CDM model in which we 
assume an equation of state evolving with redshift, the 
f^CDM model where we allow the spatial curvature of 
the universe to vary, and the model in which an interac- 
tion among the dark matter and dark energy sectors is 
switched on, the £CDM model. These scenarios are an 
extension of the minimal cosmological model plus three 
(N Us ) active (sterile) massive neutrino species. We con- 
sider subsets of the following parameters: 



{uj b , ui c , 9 S , r, n s , log[10 10 A s ], m v , m Ve , N v „ , w(w ), w a ,Q k ,€} , 



where u b = fi b /i 2 and uj c = Q c h 2 are the physical baryon 
and cold dark matter densities, S is the ratio between 
the sound horizon and the angular diameter distance at 
decoupling, r is the optical depth, n s is the scalar spectral 
index, A s is the amplitude of the primordial spectrum, 
m„ is the active neutrino mass, m Vs is the sterile neutrino 
mass and N Va is the number of thermalized sterile neu- 
trino species 2 , w is the dark energy equation of the state 
parameter, f2j. is the curvature parameter, £ is the dimen- 
sionlcss parameter which encodes the dark matter dark 
energy interaction, and wq, w a are parameters related to 
the dark energy equation of state. Table Q] specifics the 
priors considered on the different cosmological parame- 
ters. In all cases in which 111 / -1 we have considered 
the effect of dark energy perturbations, fixing the dark 
energy speed of sound c 2 = 1 . 

Our basic data set is the seven-year WMAP data 
[H [HI (temperature and polarization) with the routine 
for computing the likelihood supplied by the WMAP 
team. We consider two cases: we first analyze the 
WMAP data together with the luminous red galaxy 
clustering results from SDSS-II (Sloan Digital Sky Sur- 
vey) [32| and with a prior on the Hubble constant from 
HST (Hubble Space Telescope) [HI], referring to it as 
the "runl" case. We then include with these data sets 
Supernova la Union Compilation 2 data [34j . and we 
will refer to this case as "run2". In addition, we also 
add to the previous data sets the BBN measurements of 



the 4 He abundance, considering separately helium frac- 
tions of Fp 1 = 0.2561 ± 0.0108 (see Ref. ||) and of 
Y p 2 = 0.2565 ± 0.0010 (stat.) ±0.0050 (syst.) from 
Ref. [36[. Finally, we also consider the Deuterium abun- 
dance measurements \og(D/H) = —4.56 ± 0.04 from 
Ref. [33. 

Since the CMB data we are considering are not sig- 
nificantly constraining the amount of primordial He- 
lium abundance, we fix it to the value Y p = 0.24, 
consistent with current observations. Then we use the 
MCMC chains from each different run and perform im- 
portance sampling obtaining the predicted values for Y p 
and \og(D/H) with an interpolation routine using a gri d 
of the public available PArthENoPE BBN code (see [38j]) 
for each point (wb, N° e = 3+N Vs ) of a given cosmological 
model, as in [39j. 

In the following, we will present the cosmological con- 
straints on the masses of the active and the sterile neu- 
trino states as well as on the number of sterile states for 
different cosmological scenarios, namely, a universe with 
a constant equation of state w ^ — 1, a universe with 
a time varying dark energy fluid, a universe with a non 
vanishing curvature component and a universe with in- 
teracting dark matter-dark energy sectors. 



A. 111CDM cosmology 



2 We assume that both active and sterile neutrinos have a degen- 
erate mass spectra. 



We first consider a cosmological model including stan- 
dard cold dark matter and a dark energy fluid character- 
ized by a constant equation of state w . We consider the 
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Parameter 


Prior 




0.005 - 


-> 0.1 


tt c h 2 


0.01 - 


► 0.99 


e s 


0.5 - 


> 10 


T 


0.01 - 


-»■ 0.8 


n s 


0.5 - 


> 1.5 


ln(10 1( M s ) 


2.7 - 


-> 4 


m Vs [eV] 


- 


> 3 


m v [eV] 


- 


► 3 




- 


► 6 


w(wo) 


-2 - 


-> 


w a 


-1 - 


-> 1 




-0.02 - 


-+ 0.03 




-2 - 


-»■ 0. 



TABLE I: Flat priors for the cosmological parameters consid- 
ered here. 



following set of parameters: 

{uj b , u c , S , t, n s , log[10 10 A s ], m„, m„, , A^ s , w} 



(1) 



Table [TT| shows the one dimensional (ID) marginalized 
95% CL bounds on N Vs , m Vs and m v using the two com- 
binations of data sets described above. Note that the 
addition of SNIa data affects only the number of mas- 
sive sterile neutrino species, and not their masses. The 
bounds obtained in a ACDM scenario, see Tab. I VII in Ap- 
pendix [A] [iU are slightly relaxed when the dark energy 
equation of state is allowed to vary. There is a strong 
and very well known degeneracy in the m v — w plane 
(and therefore also in the m Va — w plane) as first no- 
ticed in Ref. (23|. Cosmological neutrino mass bounds 
become weaker if the dark energy equation of state is 
taken as a free parameter. If w is allowed to vary, the 
cold dark matter mass energy density f2 c can take very 
high values, as required when m„ (or m Ue ) is increased 
in order to have the same matter power spectrum. We 
can observe this degeneracy in the upper panel of Fig. [T] 
There exists also a degeneracy between the number of 
sterile neutrino species and the dark energy equation of 
state, see Fig. Q] (lower panel). Sub-eV massive sterile 
neutrino species may be quasi relativistic states at de- 
coupling. One of the main effects of N Vs comes from the 
change of the epoch of the radiation matter equality, and 
consequently, from the shift of the CMB acoustic peaks, 
see Ref. |4(J for a detailed study. The position of acoustic 
peaks is given by the so-called acoustic scale 8 a, which 
reads 



0a = 



(^rec) 



(2) 



where rg(z rec ) and r s (z rec ) are the comoving angular di- 
ameter distance to the last scattering surface and the 
sound horizon at the recombination epoch z rec , respec- 
tively. Although rs(z rec ) almost remains the same for 
different values of N Vs , r s (z rec ) becomes smaller when 
Af„ s is increased. Thus the positions of acoustic peaks 



are shifted to higher multipoles (smaller angular scales) 
by increasing the value of N Vg [41[. A dark energy com- 
ponent with w > — 1 will decrease the comoving angular 
diameter distance to the last scattering surface r$(z rec ), 
shifting the positions of the CMB acoustic peaks to larger 
angular scales, i.e. to lower multipoles £, compensating, 
therefore, the effect induced by an increase of N v . 




FIG. 1: The upper panel shows the 68% and 95% CL bounds 
from "run 1" in the plane w — m v (in blue) and w-m Vs (in red), 
respectively. The masses of the sterile and active neutrinos 
are both in eV. The lower panel depicts the 68% and 95% 
CL contours in the plane w — N„ s for a constant dark energy 
equation of state. 



We have also computed the constraints after combin- 
ing the results of "run 2" with those coming from differ- 
ent measurements of BBN light element abundances, see 
Tab. |TTJ The bounds on N v , obtained in a wCDM cos- 
mology when adding BBN constraints are very similar to 
those obtained in a ACDM scenario [22]. The limits on 
the active and neutrino masses are mildly relaxed. 



B. to(o)CDM cosmologies 

We also consider a time varying equation of state with 
a parameterization that has been extensively explored in 
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TABLE II: ID marginalized 95% CL bounds on N Va , m Vs and m„ using the two combinations of data sets described in the 
text (rl refers to "run 1" and r2 refers to "run 2" , respectively) for the raCDM cosmology. We also show the constraints after 
combining the results of "run 2" with those coming from different measurements of BBN light element abundances. 



Parameter 


95% CL (rl) 


95% CL (r2) 


1 V 


Y' 2 




YJ + D 


Nu a 


< 4.4 


< 3.9 


< 2.3 


< 1.6 


< 1.6 


< 1.3 


m v [eV] 


< 0.34 


< 0.31 


< 0.25 


< 0.23 


< 0.23 


< 0.23 


m„ s [eV] 


< 0.51 


< 0.57 


< 0.68 


< 0.71 


< 0.75 


< 0.75 



the literature (43-|45|: 

w(a) = w Q + w a (l - a) . (3) 

We consider the following set of parameters: 

{uj b , uj c , 6 S , t, n s , log[10 10 A,] , m v , m Us , N Vs ,w ,w a } . 

(4) 

Table Hill shows the ID marginalized 95% CL bounds on 
N„ s , m Vs and m„ using the two combinations of data sets 
used along this manuscript. The addition of SNIa data 
does not improve at all the constraints from "run 1". 
There exist large degeneracies in the m u — wq, m lJs — wo 
and N Vs — wq planes. These degeneracies are identical to 
the wCDM cosmology ones, and therefore we will not il- 
lustrate them here to avoid redundancy. We show instead 
the mild degeneracy in the N v> — w a plane, see Fig. [2] 
In this case, an increase of N Vs will be compensated by 
a decrease of w a - The BBN bounds combined with the 
"run 2" constraints on the neutrino parameters are also 
depicted in Tab. [TTTJ 
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w 

a 

FIG. 2: 68% and 95% CL contours arising from "run 1" in 
the plane w — N Vs for the ra(a)CDM cosmology. 



C. Non flat SJtCDM cosmologies 

We also explore here the constraints in non flat cos- 
mologies. We consider the following set of parameters: 

{w b ,w c ,6 s ,r, n s ,\ogll0 10 A s ],m u ,m„ 3 ,N„ 3 ,tt k } . (5) 



Current CMB measurements combined with SNIa and 
BAO data, give the constraint Q k = -0.0057l£;^ @]. 
Table HVl shows the ID marginalized 95% CL bounds on 
N Vs , m Vs and m„ using the two combinations of data sets 
considered here. The bounds are very similar to those 
obtained in the other two cosmological models explored 
above. 

Figure [3] shows the existing degeneracies between the 
curvature energy density and the number of massive ster- 
ile neutrino species N Ug and between the curvature and 
the masses of the active and sterile neutrino species. In 
an open universe with fife > the comoving angular di- 
ameter distance to the last scattering surface rg(z rec ) will 
be higher, shifting the positions of the CMB acoustic 
peaks to smaller angular scales, i.e. to larger multipoles 
£, which can be compensated with a larger dark mat- 
ter energy density which will allow for higher neutrino 
masses. At the same time a higher dark matter energy 
density will decrease the height of the acoustic peaks, 
features which can be compensated by a larger number 
of sterile neutrino species. 

Table HVl also show the constraints arising from BBN 
measurements combined with "run 2" bounds. 



D. Coupled £CDM cosmologies and modified 
gravity scenarios 

Interactions within the dark sectors, i.e. between cold 
dark matter and dark energy, are still allowed by ob- 
servations (see Ref. [46| and references therein). We 
parametrize the dark matter-dark energy interactions at 
the level of the stress-energy tensor conservation equa- 
tions, introducing an energy momentum exchange of the 
following form [47] : 

= -Qu , (6) 

with 

Q v = ZHpdeul/a or Q v = £Hpd e ut e /a , (7) 

where u^ de ^ is the cold dark matter (dark energy) four 
velocity and £ is a dimensionless coupling, considered 
negative in order to avoid early time non adiabatic in- 
stabilities [27]. In general, coupled models with Q u pro- 
portional to uf, e are effectively modified gravity models. 
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TABLE III: As Tab. HU but for the iu(a)CDM cosmology. 



Parameter 


95% CL(rl) 


95% CL(r2) 


Y l 
1 p 


Y' 2 
1 p 


y; + d 






< 4.5 


< 4.4 


< 2.6 


< 1.7 


< 1.7 


< 1.4 


m v [eV] 


< 0.33 


< 0.31 


< 0.26 


< 0.23 


< 0.22 


< 0.22 


m„ s [eV] 


< 0.49 


< 0.48 


< 0.58 


< 0.66 


< 0.73 


< 0.72 



TABLE IV: As Tab.|I[] but for the fi fe CDM cosmology. 



Parameter 


95% CL(rl) 


95% CL(r2) 


i p 


Y' 2 

p 


y; + d 


Y*+D 




< 4.3 


< 3.8 


< 2.4 


< 1.7 


< 1.7 


< 1.4 


m v [eV] 


< 0.34 


< 0.31 


< 0.24 


< 0.22 


< 0.22 


< 0.21 


m„ s [eV] 


< 0.45 


< 0.52 


< 0.70 


< 0.72 


< 0.77 


< 0.75 




-0.02 



-0.01 



0.01 0.02 




► " 0.3 



FIG. 3: The upper panel shows the 68% and 95% CL bounds 
from "run 1" in the flk — N Vs plane. The lower plane shows 
the analogous in the fi/j — m v (in blue) and fij. — m Va (in 
red) planes, respectively. The masses of the sterile and active 
neutrinos are both in eV. 



We consider the following set of parameters: 
{ub,Uc,Os,T,n s ,log[10 10 A s ],m„,m,v B ,N l/s ,w,£,} . (8) 
Interactions between the dark matter and dark en- 



ergy sectors can relax the bounds on the active neutrino 
masses [26|, [27| since, for negative couplings, the power 
spectrum increases due the the higher matter energy den- 
sity when £ < 0. Such a power enhancement effect (in- 
duced by the presence of a coupling) can be compensated 
by adding massive neutrinos in the game, which will sup- 
press the matter power spectrum. Therefore, there exists 
a well-known m v — £ degeneracy. Here we illustrate an 
additional degeneracy, the one existing in the N Vs — £ 
plane. Figure 0] shows the existing degeneracy among 
the coupling £ and the number of sterile massive neutrino 
species for a coupling term proportional to u c v (identical 
results are obtained for the case proportional to uf, e ). As 
first noticed in Ref. [27J a huge degeneracy is present be- 
tween the coupling £ and the mass energy density of cold 
dark matter u) c , with the former two quantities having 
positive correlations. In a universe with a negative dark 
coupling £, the matter content in the past is higher than 
in a standard ACDM scenario, since the cold dark matter 
and dark energy densities read 



(w + l) 



(9) 
(10) 



with ft® de the current cold dark matter (dark energy) 
mass energy densities. Therefore, the amount of intrinsic 
dark matter needed - that is, not including the contribu- 
tion of dark energy through the coupling term - should 
decrease as the dark coupling becomes more and more 
negative. Therefore, the number of effective neutrino 
species should also decrease as the coupling gets more 
negative to leave unchanged both the matter-radiation 
equality epoch and the first CMB peak heigh. Table 
IVl shows the analogous to previous sections but for the 
coupled case. Notice that the bounds on the number of 
massive sterile neutrinos arising from "runl" and "run 2" 
analyses are milder than those obtained in the previous 
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cosmologies, due to the degeneracy between the coupling 
£ and N Ve . 




FIG. 4: 68% and 95% CL bounds arising from the "run 1" 
analysis in the plane £ — iV„ s for a universe with an interacting 
dark matter-dark energy fluid. 



III. DISCUSSION AND CONCLUSIONS 

Neutrino oscillation experiments indicate that neutri- 
nos have non zero masses and open the possibility for 
a number of extra sterile neutrino species. LSND and 
MiniBooNE antineutrino data require these extra ster- 
ile species to be massive. Much effort has been devoted 
in the literature to constrain the so called (3+1) (three 
active plus one sterile) and (3+2) (three active plus two 
sterile). Recently, global fit analyses incorporating new 
reactor antineutrino fluxes have shown that (3+2) mod- 
els with 0.5 — 1 eV sterile species provide a very good fit to 
short baseline data. Cosmology can set bounds on both 
the active and sterile neutrino masses as well as on the 
number of sterile neutrino species. We have explored here 
the current constraints on these parameters in natural ex- 
tensions of the minimal ACDM cosmology Namely, we 
have explored the neutrino constraints in scenarios with- 
out a cosmological constant as the dark energy fluid, with 
a non vanishing curvature, or with coupled dark matter- 
dark energy fluids. Figure [5] summarizes our results in 
the m v — N Ug and N Vs — m Ug planes. Notice that mod- 
els with two massive 0.5 — 1 eV sterile neutrinos plus 
three sub-eV active states are perfectly allowed at 95% 
CL by current Cosmic Microwave Background, galaxy 
clustering and Supernovae la data. Interestingly, these 
models are precisely the ones which, with the new reac- 
tor fluxes prediction, improve considerably the global fit 
to short baseline data. While we have not checked di- 
rectly the results of Ref. [28| in which models with non 
zero curvature and two extra sterile neutrinos seem to 
exclude w — —1 at 95% CL, it seems plausible that non 
standard ACDM cosmologies with sterile neutrino species 
provide also a very good fit to cosmological data. Big 













// ,4 / 



0.1 



0.2 0.3 
m 



0.4 



0.5 




FIG. 5: The upper panel shows the 68% and 95% CL bounds 
from "run 1" in the m„ — N Vs plane. The lower plane shows 
the analogous in the N Vs —ra Va plane. The masses of the ster- 
ile and active neutrinos are both in eV. Red, blue, magenta 
and green contours denote mjCDM, io(a)CDM, fifeCDM and 
£CDM cosmologies, respectively. 



Bang Nucleosynthesis constraints could compromise the 
viability of these models if the two sterile neutrino states 
are fully thermal. If the BBN constraints are obtained 
using the helium fraction measurements Y p from [35j ex- 
clusively, two extra sterile neutrino states are perfectly 
allowed (N Vs < 2.3 at 95% CL). Even after combining 
these data with Deuterium measurements, two massive 
neutrino states are only marginally excluded (iV„ 3 < 1.7 
at 95% CL) . The tightest bound on the number of sterile 
neutrino species arises when helium measurements from 
Ref. [36j are combined with Deuterium data. Further de- 
velopments in BBN determinations of light element abun- 
dances may have a large impact in further constraining 
the number of sterile neutrino species N Vs . 
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TABLE V: As Tab. |H] but for the £CDM cosmology. 



Parameter 


95% CL(rl) 


95% CL(r2) 


Y l 
1 p 


Y' 2 
1 p 


YJ+D 


YJ+D 




< 5.2 


< 4.7 


< 2.4 


< 1.7 


< 1.8 


< 1.4 


m v [eV] 


< 0.34 


< 0.33 


< 0.21 


< 0.21 


< 0.20 


< 0.20 


m„ a [eV] 


< 0.35 


< 0.38 


< 0.47 


< 0.51 


< 0.49 


< 0.49 
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parison with the results obtained within the other cos- 
mological scenarios explored here. 



Appendix A: ACDM results 

Tabic I VII summarizes the main results obtained in [1H 
for a standard ACDM cosmology, for the shake of com- 
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